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Abstract 

Consider the Delaunay triangulation T of a set P of points in the plane as a Euclidean 
graph, in which the weight of every edge is its length. It has long been conjectured that the 
dilation in T of any pair p,p' G P, which is the ratio of the length of the shortest path from 
p to p' in T over the Euclidean distance can be at most tt/2 rss 1.5708. In this paper, 

we show how to construct point sets in convex position with dilation > 1.5810 and in general 
,— i position with dilation > 1.5846. Furthermore, we show that a sufficiently large set of points 

drawn independently from any distribution will in the limit approach the worst-case dilation for 
that distribution. 



C/2 



1 Introduction 



In this paper we establish new lower bounds for the dilation of the Delaunay triangulation of points 

ON 

in the plane. We say that a graph G = (P, E) on a set P of points in the plane is a Euclidean graph 
\q if the weight of each edge (p, q) € E is the Euclidean distance \\pq\\. The dilation for any pair of 

points p,p' G P is the ratio of the length of a shortest path from p to p' in G over the Euclidean 



distance The dilation of G is the maximum dilation between any pair. The dilation of G 

has also been called the stretch factor or spanning ratio of T, and the concept has been used to 



> 

analyze routing algorithms for networks [2, 10, 17] and to define spanners [8, 9, 13-15]: a t-spanner 
is a graph defined on a set of points such that the dilation between any two points is at most t. 

A Delaunay triangulation of P is a triangulation of the convex hull of P in which the circumcircle 
of every triangle contains no points of P in its interior. Note that each circumcircle will have at 
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Figure 1: The dilation of this Delaunay triangulation of a set of points on a circle approaches tt/2. 

least the three vertices of the triangle on its boundary, and that if it has more points of P on the 
boundary, the Delaunay triangulation is not unique. Our definition allows any triangulation of the 
points on a circle with empty interior to be called a Delaunay triangulation. 

One of the first results in computational geometry on spanners was a proof of Chew's conjec- 
ture [3] that the Delaunay triangulation is a spanner in the plane. Dobkin, Friedman and Supowit 
proved that the dilation of the Delaunay triangulation of any set of points in the plane is at most 
(1 + \/5)tt/2 < 5.084 [7]; the best upper bound known is t = (4^/3/9)vr < 2.419 by Keil and 
Gutwin [12]. At CCCG 2007 the first author posed, as an open problem [6], improving dilation 
bounds for the special case of a set of points in convex position. In this issue, Cui, Kanj and Xia 
establish an upper bound of 2.33 for points in convex position [5], and we establish a new lower 
bound. 

Until now, the only lower bound on the dilation for Delaunay triangulations was Chew's con- 
struction, which achieves tt/2 — e, for any e, by sampling points P uniformly on a unit circle [3]. 
Identify two antipodal points p and p', and triangulate P by taking all edges nearly perpendicular 
to pp', as in Figure 1. Since all points are co-circular, any planar triangulation of P is a valid 
Delaunay triangulation; alternatively, we could perturb the points to break the co-circularity, and 
make this the only triangulation. The shortest path from p to p' via the triangulation follows the 
boundary of the circle, so its length approaches ir. On the other hand, the Euclidean distance is 
clearly 2, so the dilation t —> tt/2. 

In their book on spanners, Narasimhan and Smid mention that "it is widely believed that, for 
every set of points in R 2 , the Delaunay triangulation is a (7r/2)-spanner" [15]. We show here how 
to construct point sets in convex position whose Delaunay triangulation is not a (-7r/2)-spanner, but 
in fact has dilation > 1.5810. We then modify the construction to create point sets not in convex 



2 



Figure 2: Observation 1 compares two ways to go from p to q: direct by arc of length j3 or via q' 
by arc+segment of length 9 — (3 + 2sin(0/2). These balance if = 9/2 + sin(9/2). When qq' is a 
diameter, /3 = \ + ir/2 radians. 

position whose Delaunay triangulation has slightly larger dilation; we do not yet know whether the 
maximum dilation for Delaunay triangulations can be approached by points in convex position. 
Finally, based on these constructions, we prove that as you draw sufficiently many random points 
from any single distribution, the dilation of the Delaunay triangulation approaches the worst-case 
dilation for that distribution. 

2 Constructions with dilation greater than tt/2 

To construct examples of points with dilation greater than tt/2, we use a simple observation about 
paths around sectors of unit circles that is illustrated in Figure 2: 

Observation 1 Consider the segment of a unit circle defined by points q, q' that subtend angle 6, 
and place p between them so that arc pq has angle /3. Let A be the path consisting of the arc from p 
to q' followed by the segment qq' . The arc from p to q around the boundary of this sector is shorter 
than the length of A when (3 < 9/2 + sin(0/2). 

PROOF. The arc pq has length fj, while the arc pq' and the segment qq' together have length 
9 - (3 + 2sin(#/2). 

We are now ready to describe the construction for points in convex position. Form a convex 
region bounded by two unit semicircles having centers on the x axis separated by distance d. As in 
Figure 3(a), introduce points on these semicircles uniformly and identify two points p and p' at an 
angle of a from the x-axis. Next, triangulate the semicircle with p by adding chords to the convex 
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Figure 3: (a) The basic construction consists of evenly sampled points on two unit semicircles with 
centers on the x axis separated by d. We mark two points p and p' that make an angle of a with 
the x-axis. (b) We choose the Delaunay triangulation in order to maximize the shortest path in the 
triangulation from p to p' . (c) One locally shortest path (solid red) simply follows the boundary 
and has length approaching ir + d; another (dashed) crosses the construction once and has length 
approaching ir + 2 — 2a + d. 

hull in a way that ensures that any shortest path from p to the endpoints of the semicircle follows 
the boundary of the circle. One possibility is shown in Figure 3(b). Observation 1 ensures that the 
arc is the shortest path from p to either endpoint of each chord used to triangulate. 

In our triangulation there are two types of locally optimal paths from p to p' , drawn solid 
and dashed in Figure 3(c). The first type, which follows the perimeter of the region (clockwise or 
counter-clockwise), has length w + d, since we walk around one semicircle and bridge the gap of 
width d. The other type, which crosses over via one of the vertical edges, has length 2 • (vr/2 — a) 
for the two circular arcs plus 2 + d for the straight parts, so it + 2 — 2a + d in total. Observation 1 
ensures that any other path will be longer. If we set a = 1 radian, these two path types have equal 
lengths. 

Finally, we have to compute the Euclidean distance between p and p', which is i = y/4 + d 2 + 4dcos 1. 
The dilation approaches t = (ir + d)/£. As a function of d, t remains above 1.5810528 for 
d G [0.293,0.294]. For d = 0.29, we have t > 1.581 > f . 

Theorem 2 There exists a set P of points in convex position in the plane, such that the Delaunay 
triangulation of P has a dilation of at least 1.5810. 

In this construction 18 points suffice to give a dilation > tt/2 and 222 points give a dilation > 1.5810. 
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Figure 4: Increasing the dilation by placing points on three circles, with four shield points. Two 
unit radius circles are separated by d = 0.58, and the third circle, C, of radius 1.15057 is centered at 
their midpoint. Points p and p' , placed by Observation 1, have distance 2.4, and achieve a dilation 
of 1.5846. 

If we allow points that are not in convex position, we can modify this construction, as illustrated 
in Figure 4, and increase the dilation slightly. The idea is to replace the two straight segments with 
points that lie on a common circle, C. The new triangulation edges inside the polygon will not 
be used by any shortest path from p to p'; we need to prevent short-cuts outside the polygon by 
the edges needed to complete the triangulation there. At each of the four locations where a unit 
circle meets C, we add a shield point s on the line from the center of the unit circle through the 
intersection, and draw a ray from s through the center of C. We place points densely on the arcs 
of all three circles, leaving gaps in the angles formed by the rays from each shield point s. The 
position of s on its line is chosen so that the tangents from s to the two circles form a path that is 
just longer than the path that follows the circle boundaries. The triangulation outside the circles 
is completed by fans from each shield point s. 

For the best ratio, we separate the unit circles by d = 0.58, and place circle C of radius 
r = 1.1507 midway between their centers. As depicted in Figure 4, the arcs of the unit circles are 
no longer semicircles, but subtend 29 = 2.2895 radians; the two arcs of the larger circle C each 
subtend 2/3 = 1.30432 radians, with a shielded gap of g = 0.0065. When we place p and p' by 
Observation 1 we obtain a straight line of length I = 2.4 and a dilation (8 + r/3 + g)/£ > 1.5846. 

Theorem 3 There exists a set P of points in the plane such that the Delaunay triangulation of P 
has a dilation of 1.5846. 
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3 The dilation of the Delaunay for random point sets 

Because a small number of points may determine the dilation of a large Delaunay triangulation, 
we observe that dilation factors greater than ir/2 are actually the rule, and not the exception, 
for sufficiently large sets of random points. In this section we formally prove this for quite weak 
restrictions on the probability density function. 

Although the constructions of the previous section create degenerate point sets, it is not hard 
to alter them slightly, for any finite number of points k, so that the points are in general position 
and the dilation is nearly the same. Once the points are in general position, there is a value 5 > 
such that perturbing every point by at most S leaves the combinatorial structure of the Delaunay 
triangulation unchanged. (Abellenas et al. [1] provide an algorithm to compute 5.) The idea of this 
section is to show that whenever there is some finite probability that k points chosen at random 
will create a Delaunay configuration with dilation greater than ir/2, then when we choose n 3> k 
points, we will, with high probability, find such a /c-configuration. 

The material below applies to general settings in which we seek a lower bound L*, that applies 
with high probability as n — > oo, on a function L n on n points in R d that is a maximum over 
configurations in products of intervals that are subsets of R d . Assume the following two conditions: 

(i) For each n, L n is scale and translation invariant, that is, L n (axi+b, . . . , ax n +b) = L n (x±, . . . , x n ) 
for all a ^ 0, b € M d . 

(ii) There exists a value L* such that for all e > 0, there exists a 5 > 0, integer K > 1, and points 
x\, . . . , xk € [0, l] d , such that if B(x, 6) denotes the ball of radius 5 in we have 

inf inf inf L K+m (z 1 , . . . , z K ,yx, . . . ,y m ) > L* - e. (1) 

zieB( X1 ,S),...,ZKeB{xK,S) m yi,...,y m £l0,l] d 

By way of remark, the easiest way to satisfy (ii) is for L n to be a continuous function of its inputs; 
condition (1) would then follow if the value L* were achieved for a certain point set. Unfortunately, 
the dilation of the Delaunay triangulation is not continuous; moving points by a small amount may 
induce a diagonal flip that changes the shortest path length. Thanks to Abellenas, Hurtado, and 
Ramos [1], however, we do have (ii) for non-degenerate point sets. Briefly, (ii) is satisfied if there 
exists a small nonempty ball in the neighborhood of (but not necessarily covering) a point at which 
L* is reached on which L n is bounded from below by L* — e. 
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To conclude this paper, we prove a general theorem on functions L n that satisfy our two 
conditions. It essentially says that a copy of some pessimal construction can be expected in a 
random point set under weak restrictions. The probability density / mentioned in the theorem 
is arbitrary (i.e., is the derivative of an arbitrary continuous cumulative distribution function). It 
can have unbounded support and infinite peaks, and could possibly fail to be continuous at almost 
all x. 

Theorem 4 Let L n : (TR d ) n — > IR satisfy (i) and (ii) for a certain value L* . Let Xi,...,X n be 
i.i.d. random vectors drawn from a common density f in TR d , then for every e > 0, 

lim V{L n (X 1 ,...,X n ) > L* — e} = 1. 

3.1 Preliminaries 

Let S C [0, l] d be a fixed set of positive volume. Let / be a density of IR d . We say that x is a 
Lebesgue point for S if the density does not vary wildly in the neighborhood of x, expressed as: 

lim ^rsf(y)dy = lim Usf(v)*v = m 

r +° Jx+rS d y ri0 r d J s dy 
Thus, for a Lebesgue point x, the average density value over a small ball centered at x tends to 
the value of the density at x as the radius decreases. A fundamental property of densities is that 
almost all points are Lebesgue points. 

Lemma 5 For any density f , and any set S C [0, l] d of positive volume, almost all x are Lebesgue 
points for S. 

Proof. See, e.g., Wheeden and Zygmund [16], p. 108, or de Guzman [11]. □ 

Unfortunately, we cannot conclude that the density is near uniform around Lebesgue points; 
convergence of averages does not imply uniform convergence. In particular, it could be that at all 
x, the infimum of the density over any such ball of any positive radius is zero, and the supremum 
is infinite. Therefore, additional technical work is needed. We extend the notion to sets of points 
x that are Lebesgue points for a finite collection S\, . . . , Sk from [0, l] d . Clearly, almost all x are 
still in this set. The next Lemma says that we can always find many Lebesgue points that are far 
apart and have certain restrictions on the densities. 
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Lemma 6 For any density f, and any sets Si,...,Sk Q [0,1] of positive volume, there exist 
constants < a < b < oo such that for any integer N > 0, we can find So > 0, xi, . . . , xjy £ IR d , 
and constant c > suc/i i/tai 6 > /(x«) > a, i/ie distance between each Xi and xj, is more than 
c, and 

for all I <i < N, < j < K, < S < S , where S = [0, V\ d . 

PROOF. The existence of the claimed point set is shown by the first moment method. Consider 
i.i.d. random points X[, X' 2 , . . . drawn from /. Choose a and b such that V{f(X[) [a, b]} < 1/2. 
Let Xi, . . . , Xjy be the first N of the X'j's that have f(Xi) 6 [a, b\. Let T be the number of random 
vectors needed to achieve this. Note that £{T} < 2N, as T is stochastically smaller than a sum 
of n geometric random variables with parameter 1/2. With probability one, all X'^s, and thus all 
X,'s, are Lebesgue points for So = [0, l] d , Si,..., Sk- Thus, there exists So (depending upon N, f, 
Xi,..., X N ) such that for S < S , for all 1 < i < N, < j < K, 

where A denotes Lebesgue measure. The probability that the distance between Xi and Xj is 
less than c for some i ^ j is not larger than £{T 2 M(c)}, by the union bound, where 

M(t) ^ sup / f(y) dy. 

x Jx+t[0,l] d 

As £{T 2 } = 0(N 2 ), this is less than 1/2 by choice of c > because M(c) I as c I 0. We take such 
a c (depending upon N). Thus, by the first moment method, there exists a So,x\, . . . ,xjy,c > 
such that jointly b > f(xi) > a, the distance between each x% and Xj is more than c, and 

f4>L s / €[/(x,)/2 ' 2/(x,)1, 

for all 1 < i < N, < j < K, < 5 < S - □ 
3.2 Proof of Theorem 4 

We are now ready for the proof of Theorem 4, which proceeds by the second moment method. Let 
e > be given. Find a finite K = K e , points xi, . . . ,xk in [1/4, 3/4] d and r = r e > such that 
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(1) holds for any vector . . . , x' K ) with x\ G B(x{,r), 1 < i < K . In our application, this is a 
locally stable Delaunay configuration with large dilation. Define Si = B(xi,r). By choice of r, we 
can insure that the Si are non-overlapping and completely contained in [0, l] d . 

Fix a and b as in Lemma 6. Armed with the sets Si, and given any integer N > 0, Lemma 6 
tells us that we can find 5q > 0, yi, . . . , yx £ IR d , and constant c > such that b > f(xi) > a, the 
Lqo distance between each Xi and Xj is more than c, and 

for all 1 < i < N, < j < K, < 6 < S , where S = [0, l] d . 

Choose 5 = l/n 1 /^, and let n be large enough to insure that 5 < min(c, 5o). For a set A, we 
define its cardinality by the counting measure, 

n 

\ A \=J2 1 [x i eA], 

i=i 

where X±, . . . , X n are the i.i.d. data drawn from /. For 1 < i < N, define the events 

Ei = nf =1 [\ Xi + SSj\ = i] n [\ Xi + S([o, if - u^kSj)] = o 

It is clear that if one of the events Ei happens, then 

L n (Xi, . . . ,X n ) > L* - e. 
Set pi = V{Ei} and pij = V{Ei n Ej}. We have, by the second moment method [4], 

V{UiEi}> = ( ^ )2 . 

It helps to consider the ratios 

Pij 
PiPj ' 

Define 

qi{£) = f f,l<£<K, 
J Xi+6Si 

and qi(0) = / x . +< j[ i]d /• Then by the multinomial formula, 

n - 1: v . . n 1Xn _ K y.W ■ ■ ■ «(i0 (i - »(0))-» = ^ »(i) • ■ ■ (i - «(0))-» 
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and 



Thus, 



Vi3 = (n-2K)\ qi{1) " " " * {KM1) " " " V {K) (1 " * (0) " q ^~ 2K 



Pij (n- K)\{n - K)\ (1 - 9i (0) - qMT~ 2K 



PiPj n\{n-2K)\ (i _ qi (0)) n - K (1 - 9j (0)) n " A 
(n-K)!(n-K)! . . , . VX _ K 



1 - 46J d ) K 

-K 



= (1 - 46/n) 
1 

< 



_ AbK 

dcf 



1 - 

n 
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if n > 46.fr. Thus, we have 

> 



2 



CLiPi) 2 

J2iPi + i(J2iPi) 2 
i 

7 + iHiPiY 1 

We have 7 — >■ 1 as n — > 00. Furthermore, can be made as large as desired. To see this, note 

that ft (0) < 2f{ Xi )5 d < 2b/n, qi{i) > (1/2) f(xi)5 d \{St) > (l/2)a\(S e )/n d = £/n for t > as each 
has the same positive volume. Thus, 

2> 2 A '(^W« / " )t ' (1 " 2hln) "~" ~ ~^T~~' 

i 

which is as large as desired by choice of iV. Therefore, 

lim T{L n {X u . . . ,X n ) > L* - e} = 1. □ 
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4 Open problems 



There remains a large gap between the upper and lower bounds on the dilation of the Delaunay 
triangulation, both for points in convex and general position. We still believe that the true bound 
is near 1.6 in both cases; even in the convex case it would be interesting to reduce the upper bound. 
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